Nonpersistence of resonant caustics in perturbed elliptic 
billiards 



^ ■ Sonia Pinto-de-Carvalhof and Rafael Ranurez-Rosf 

f Departamento de Matematica, ICEx, Universidade Federal de Minas Gerais, CP 702, 
30123-970, Belo Horizonte, MG, Brazil 

| Departament de Matematica Aplicada I, Universitat Politecnica de Catalunya, Diagonal 
'■ 647, 08028 Barcelona, Spain 

(N 

Abstract. Caustics are curves with the property that a billiard trajectory, once tangent to it, 
' stays tangent after every reflection at the boundary of the billiard table. When the billiard table 

£^ . is an ellipse, any nonsingular billiard trajectory has a caustic, which can be either a confocal 

ellipse or a confocal hyperbola. Resonant caustics — the ones whose tangent trajectories are 
closed polygons — are destroyed under generic perturbations of the billiard table. We prove 
that none of the resonant elliptical caustics persists under a large class of explicit perturbations 
of the original ellipse. This result follows from a standard Melnikov argument and the analysis 
of the complex singularities of certain elliptic functions. 
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1. Introduction and main result 

Birkhoff [3] introduced the problem of convex billiard tables more than 80 years ago as a way 
to describe the motion of a free particle inside a closed convex smooth curve. The particle is 
reflected at the boundary according to the law "angle of incidence equals angle of reflection". 
Good modern starting points in the literature of the billiard problem are lfTTl[T8l . 

Caustics — curves with the property that a billiard trajectory, once tangent to it, stays 
tangent after every reflection — are the most distinctive geometric objects inside billiard 
tables, since they are a geometric manifestation of the regularity of their tangent trajectories. 
For example, integrable billiards have a continuum of caustics, whereas the nonexistence 
of caustics inside a convex billiard table implies that there are some billiard trajectories 
whose past and future behaviours differ dramatically. See, for instance, [fT3l . Hence, the 
existence and persistence of caustics are two fundamental questions in billiards. Most of the 
literature deals with convex caustics, since they are easier to understand and related to ordered 
trajectories. Two exceptions are [|8l §3] and ifiQl . 

We summarize the classical existence results as follows. On the one hand, if the boundary 
curve is smooth enough and strictly convex, then there exists a collection of smooth convex 
caustics close to the boundary of the table whose union has positive area [0[T2]|. On the other 
hand, Mather [fT3l proved that there are no smooth convex caustics inside a convex billiard 
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table when its boundary curve has some flat point. Gutkin and Katok [[8]| gave a quantitative 
version of Mather's theorem. 

The robustness of a smooth convex caustic is closely related to the arithmetic properties 
of its rotation number, which measures the number of turns around the caustic per impact. 
Caustics with Diophantine rotation numbers persist under small perturbations of the boundary 
curve. This follows from standard KAM arguments 0[T2]|. On the contrary, resonant caustics 
— the ones whose tangent trajectories are closed polygons, so that their rotation numbers are 
rational — are fragile structures that generically break up. See, for instance, ifToll . 

This raises two complementary questions. First, to characterize the perturbations that 
preserve/destroy a given resonant caustic of a billiard table. Second, to determine all resonant 
caustics that are preserved/destroyed under a given perturbation of an integrable billiard table. 
These questions have been studied by several authors. Baryshnikov and Zharnitsky proved 
that the perturbations preserving a given resonant caustic of a smooth convex billiard table 
form an infinite-dimensional Hilbert manifold. As a sample, we point out that this Hilbert 
manifold is given by the set of billiard tables with constant width when the rotation number of 
the unperturbed caustic is one half ifTOl . Concerning the second question, Ramfrez-Ros lfT6l 
gave a sufficient condition for the break-up of the resonant circular caustics inside a circular 
billiard table, in terms of the Fourier coefficients of the perturbation, see Remark |3] below. 

In this paper we tackle the second question when the billiard boundary is an ellipse. In 
that case, the billiard dynamics is integrable and any billiard trajectory has a caustic [fT8l . The 
caustics are the conies confocal to the original ellipse: confocal ellipses, confocal hyperbolas, 
and the foci. Poncelet [15J showed that if a billiard trajectory inside an ellipse is a closed 
polygon, then all the billiard trajectories sharing its caustic are also closed polygons. Even 
more, if a billiard trajectory tangent to one of the elliptical caustics is a (m, n)-gon — a closed 
polygon with n sides that makes m turns around its caustic — , then all the billiard trajectories 
sharing its caustic are also (m, n)-gons, and their caustic is called (m,n)-resonant. (These 
two definitions are not restricted to billiards inside ellipses.) We shall see in Section |4] that 
there is a unique (m, n) -resonant elliptical caustic for any relatively prime integers m and n 
such that 1 < m < n/2. Our main result is that all these resonant elliptical caustics break up 
under a large class of explicit perturbations of the original ellipse, see Theorem[TJ 

The following notations are required to state the main result. Once fixed the ellipse 

f 2 2 \ 

Q = Ux,y) E R 2 : ^ + y = 1 j , a > b > 0, 

we consider its associated elliptic coordinates (//, if) given by the relations 

x = c cosh /i cos v?, y = csinh/isiny?, 

where c = y/ a 2 — b 2 is the semifocal distance of Q. The equation of the ellipse Q in this 
elliptic coordinates is fj, = fj, , where cosh /x = a/c and sinh fi = b/c. Hence, any smooth 
perturbation Q e of the ellipse Q can be written in elliptic coordinates as 

A* = Veiv) = Ho + epi(<p) + 0(e 2 ), (1) 
for some 27r-periodic smooth function fi e ((p). 
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Theorem 1. Let Hi(f) be a 2-a-periodic entire function. Iffi^ip) is not constant (respectively, 
^'i(<p) is not ix-antiperiodic), then none of the (m, n) -resonant elliptical caustics with odd n 
(respectively, even n) persists under the perturbation ([7]). 

Our proof is based on the study of the persistence of the resonant rotational invariant 
circles (resonant RICs) of some twist maps by means of a first-order Melnikov method. Only 
convex caustics can be related to the RICs of those twist maps. Thus, there is no direct way to 
extend the same procedure to the nonconvex caustic hyperbolas, but we believe that the same 
results hold for them. 

Remark 1. If fi e (ip) is constant, then the perturbed curves Q e are ellipses, so all caustics 
(resonant or not) are preserved. Hence, the hypothesis fj,i((p) nonconstant is natural, since 
we are using a first-order method. Nevertheless, we can still state some results when this 
hypothesis fails. More precisely, let us assume that 

fi e ((p) = no + e/ii -\ h e^Vi-i + eVi(^) + 0(e i+1 ), 

for some /x , • • • , IM-i £ ^ an d some nonconstant 27r-periodic entire function fJLi((p). Then: 

• If n is odd, all the (m, n)-resonant elliptical caustics with odd n break up. This result is a 
corollary of Theorem [IJ It suffices to consider 5 = e % as the new perturbative parameter, 
Q* = {fi = Ho + ■ ■ ■ + e l_1 /ij„i} as the unperturbed ellipse, and to realize that Q e is a 
0(5) -perturbation of Q* whose first-order term in 5 verifies the hypotheses of TheoremQ] 

• If n is even, we believe that all (m, n) -resonant elliptical caustics also break up, even if 
Hii^f) is 7r-antiperiodic, but we should use a second-order Melnikov method in order to 
prove it. Unfortunately, the computations become too cumbersome. 

Remark 2. If we write the perturbed ellipse Q £ in Cartesian coordinates as 

x 2 /a 2 + y 2 /b 2 + eP 1 (x,y) + 0(e 2 ) = 1, 

then 2(a 2 sin 2 y? + b 2 cos 2 tp)p,i(<p) + abP±(a cos if, b simp) = 0. In particular, the function 
Hi (</?) is 7r-antiperiodic when P 1 (x, y) is odd. 

Remark 3. The case of perturbed circular tables was studied using similar techniques in lfT6ll . 
but the final result was quite different. Let us recall it for comparison. Any billiard trajectory 
inside a circle of radius r has some concentric circle of radius ^/r 2 — X 2 as caustic, where 
< A < r plays the role of a caustic parameter. If A = r sin(m7r/n), then the circular 
caustic is (m, n)-resonant. Let us write the perturbed circle in polar coordinates (r, 9) as 

r = r € (6) = ro(l + e ri (6) + 0(6 2 )), (2) 

for some smooth function r e : T — > M. Let X^ez be the Fourier expansion of r\(6) 
and n > 2. If there exists some / G nZ \ {0} such that f\ ^ 0, then the (m, n)-resonant 
circular caustics do not persist, see ffT6l Theorem 1]. In particular, it is not known if the 
(m, ?2)-resonant circular caustics with odd (respectively, even) n break up when r\{6) is not 
constant (respectively, r[(9) is not 7r-antiperiodic). 

We complete this introduction with a note on the organization. In Section[2]we develop a 
general Melnikov theory to study the persistence of resonant RICs of twist maps. The general 
setup is adapted to billiard maps in Section Finally, Theorem \T\ is proved in Section @]by 
analysing the complex singularities of certain elliptic functions, an idea borrowed from [0. 
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2. Break-up of resonant invariant curves in twist maps 

This section is a generalization of [fT6l §2], although several hypotheses have been weakened. 
Namely, the unperturbed map can be nonintegrable, the resonant invariant circle does not 
need to be horizontal, and the shift on the invariant circles can be nonconstant. In spite of 
it, the essential idea does not change. A similar theory is contained in ifTTll . For a general 
background on twist maps we refer to the book BH §9.3] or to the review 031. 

Let T = R/27rZ, and 7Ti : T x R — > T be the natural projection. Sometimes it is 
convenient to work in the universal cover R of T. We will use the coordinates (x, y) for both 
T x R and R 2 . The lines of the form x = constant and y = constant will be called vertical 
and horizontal, respectively. A tilde will always denote the lift of a function or set to the 
universal cover. If g is a real- valued function, dig denotes the derivative with respect to the ith 
variable. We will assume that all the considered objects are smooth. Here, smooth means C°°. 
In particular, all the dependences on the perturbative parameter e are assumed to be smooth. 

We will consider certain diffeomorphisms defined on an open cylinder of the form 
Z = T x Y, for some open bounded interval Y = (y~,y+) C R. Then Z = R x Y is 
an open strip of the plane. A diffeomorphism / : Z — > Z is called an area-preserving twist 
map when it preserves area, orientation, and verifies the twist condition 

d 2 nj(x,y)^0, V(x,y)eZ. 

If the twist is positive (respectively, negative), then the first iterate of any vertical line tilts 
to the right (respectively, left). We also assume, although it is not essential, that / verifies 
some rigid boundary conditions. To be more precise, we suppose that the twist map / 
can be extended continuously to the closed cylinder T x [y_, y + ] as a rigid rotation on the 
boundaries. That is, there exist some boundary frequencies u)± E R, U- < to + , such that 
f(x,y±) = (x + w±,y±). 

Let D — {(x, x') E R 2 : w_ < x' — x < co + }. Then there exists a function h : D — > R 
such that f(x, y) = (x f , y') if and only if 

y = -dxh{x,x'), y' = d 2 h(x, x'). (3) 

The function h is called the generating function of /. Besides, if (x", y") = f(x', y'), then 

d 2 h(x,x') + d 1 h(x',x") = 0. (4) 

We study the dynamics of /, but it is often more convenient to work with the lift /, so 
we will pass between the two without comment and, in what follows, the lift / remains fixed. 

A closed curve T C Z is said to be a rotational invariant circle (RIC) of / when it is 
homotopically nontrivial and /(T) = T. Birkhoff proved that all RICs are graphs of Lipschitz 
functions. See, for instance, |fl4l §IV.C]. Let v : T — > Y be the Lipschitz function such that 
T = graphs := {(x,v(x)) : x E T}. If v is smooth, we say that T is a smooth RIC. 

Twist maps do not form a closed set under composition. For instance, the square of a 
twist map is not necessarily a twist map, and indeed typically it is not. Nevertheless, any 
power of a twist map is locally twist on its smooth RICs. 
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Lemma 2.IfT = graph v is a smooth RIC of an area-preserving twist map f : Z — >■ Z, then 

d 2 nj n (x, v{x)) ^ 0, VxeR, Vn > 1. 

iVoo/ Given any point p = (x,v(x)) G Y, let pj = (xj,v(xj)) = f j (p),tj = (l,v'(xj)), 
and Vj = (0, 1). We identify the tangent planes T P Z with the Euclidean plane R 2 . Thus, 
the vector tj is tangent to T at the point pj and Vj is a vertical vector at pj. The linear map 
df n (p) : T p Z — > T Pn Z is the composition of the linear maps df(pj) : T Pj Z — > T Pj+1 Z for 
j = 0, . . . ,n — 1. Let aj, 6j, Cj, eZ,-, a n , /3 n , 7 n , 5 n G R be the coefficients such that 

df n {p) ■ t ^ a n t n + j n v n , v ^ f3 n t n + 5 n v n . 

We note that bj = d 2 TTif(pj) and f3 n = 9 2 7Ti/ n (p). Let us suppose that the twist is positive, 
so bj > 0. We want to prove that j3 n > for any integer n > 1. The case of negative twist is 
completely analogous. 

We deduce that Cj = from the invariance of T. Hence, (3 n = ^T^j^j+h wnere 

D\ = Yll=i dk and A\ = Yii=i a k- Besides, we note that dj > because the two components 
of C \ T are invariant. Finally, we get that aj > from the preservation of orientation. □ 

Roughly speaking, a RIC is said to be resonant when all its points are periodic, but we 
need to be more precise. Let (x, y) G Z be a periodic point of the twist map /, and let n be its 
least period. Then the exists an integer m such that its lift verifies f n (x, y) = (x + 27rm, y). 
Obviously, w_ < 2nm/n < u + . Such a periodic point is said to be of type (m, n). A RIC is 
said to be (m, n)-resonant when all its points are periodic of type (m, n). 

Let / be an area-preserving twist map with a (m, n)-resonant smooth RIC T = graph v. 
Considering area-preserving twist perturbations of the form f e = f + 0(e), we prove in the 
following lemma that there exists two graphs T e = graph v t and T* = graphs* 0(e)-close 
to T and such that / e n projects the first graph onto the second one along the vertical direction. 

Lemma 3. There exist two smooth functions v e , v* : T — > Y defined for e G (— e , e ), e > 0, 
such that: 

(i) v e (x) — v(x) + O(e) and v*(x) — v(x) + O(e), uniformly in x G T; and 

(ii) f™(x,v e (x)) = (x,v*(x)),for all x G T. 

Proof. We work with the lift of the maps. Once fixed an angle x G R, let y = v(x) and 

G(y, e) := ftif?(x, y) - x - 2itm. 

This function G(y, e) verifies the hypotheses of the Implicit Function Theorem at the point 
(y,e) = (yo,0), since G(y ,0) = and d 1 G(y ,0) = d 2 7r 1 f n (x,v(x)) ^ 0, see Lemma[2l 
Consequently, there exist e , rj > such that the equation G(y, e) = has exactly one solution 
Ve = Vo + O(e) in the interval (y — r],y + rj) for all e G (— e , e ). We recall that G(y, e) had 
x G R as an extra parameter, but it appeared in a 2n -periodic smooth way. Hence, e and r] can 
be taken independent from x, the estimate \y e —yo\ = O(e) is uniform in x, and y e depends in a 
2n -periodic smooth way on x. Finally, set v e (x) = y e and then v*(x) is determined by means 
of relation f™(x, v e (x)) — (x + 27rm, v*(x)). The functions v e , v* : R — > Y are 27r-periodic 
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and smooth, so they can be projected to two smooth functions v e , v* : T — » Y that verify the 
two claimed properties by construction. □ 
We say that a (m, n)-resonant smooth RIC T of a twist map / persists under an area- 
preserving twist perturbation f t = f+0(e) whenever the perturbed map has a (m, n)-resonant 
RIC T e for any small enough e such that T e = T + 0(e). The corollary below follows 
immediately from this definition. 

Corollary 4. The resonant RIC T persists under the perturbation f e if and only ifT e = T*. 

Therefore, it is rather useful to quantify the separation between the graphs T e and T*. 
Lemma 5. v*(x) — v e (x) = L' t (x), where L e : T — >■ R is a function whose lift is 

n-1 

L e (x) = ^2h € (xj(x;e),x j+1 (x;e)), Xj(x; e) = (x,v e (x)), (5) 

3=0 

and h e is the generating function of f e . 

Proof. As long as confusion is avoided, we will omit the dependence on x and e. We 
introduce the notations (xj,§j) = fi(x,v(x)) and Wj = dxj/dx for j = 0,...,n. Then 
x = x and x n = x + 2nm, so Wq — w n — 1. Besides, y = v(x) and y n = v*(x). 
From the implicit equations ©, we get that dih(x , x\) = —y , <9 2 /i(x n _i, x n ) = y n , and 
d 2 h(xj-i, xj) + dih(xj,Xj + i) — for j = 1, . . . , n — 1. Therefore, L'(x) = dih(x , xi)w + 
Y^lZi {d2h(xj-ijXj) + dih(xj,Xj + i))wj + d2h(x n -i,x n )iD n = v*(x) — v(x). It is immediate 
to check that L : R — > K is 2it -periodic, so it can be projected to a function L : T — > R. □ 

Corollary 6. The resonant RIC T persists under the perturbation f e if and only ifL' f (x) = 0. 

We shall say that L t : T — > R is the subharmonic potential of the resonant RIC T under 
the twist perturbation f e . It is rather natural to extract information from the low-order terms 
of its expansion L € (x) = L (x) + eL\ (x) + 0(e 2 ). This is the main idea behind any Melnikov 
approach to a perturbative problem. The zero-order term L (x) is constant (and so useless), 
since L' (x) = Vq(x) — v (x) = v(x) — v(x) = 0. We shall say that the first-order term Li(x) 
is the subharmonic Melnikov potential of the resonant RIC T under the twist perturbation f e . 
The proposition below provides a closed formula for its computation. 

Proposition 7. Ifh e = h + eh 1 + 0(e 2 ), then the lift of L 1 (x) is 

n-1 

Li{x) = ^htjxj^j+t), xj = TXif 3 {x,v{x)). 

3=0 

Proof. Given any x G R, we set Xj = Xj(x) := 0) and Zj = Zj(x) := d 2 Xj(x; 0) for 
j — 0, . . . , n. Then the 0(e)-term of © is 

n— 1 

Li(x) = dih(x ,xi)z + ^2 (dih(xj,x j+1 ) + d 2 h(x j _ 1 ,x j )^z j + d 2 h(x n _ 1 ,x n )z n + 

n-1 

J2hi{xj,x j+ i). 

3=0 
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Figure 1. The billiard map f(<p, i?) = (tp', •&'). 

Using the implicit equations © for the unperturbed twist map, the first summation vanishes. 
The terms dih(x , xi)z and d 2 h(x n -i, x n )z n also vanish, since x (x; e) = x and x n (x; e) = 
x + 2iim for all e G (— e , e ). Besides, Xj = Xj(x) = Xj(x; 0) = itiP(x, v(x)). □ 
The following corollary displays the most important property of the subharmonic 
Melnikov potential in relation with the goals of this paper. 

Corollary 8. If Li{x) is not constant, then the resonant RIC T does not persist under the 
perturbation f e . 

Proof. It follows directly from Corollary [6] and the estimate L e = constant + eL! + 0(e 2 ). □ 



3. Break-up of resonant caustics in perturbed billiard tables 

Let Q be a closed strictly convex smooth curve in the plane. Let 7 : T — > Q be a 
counterclockwise parametrization. Let Z = T x (0, tt) be an open cylinder. We can model 
the billiard dynamics inside Q by means of a map / : Z — > Z, f(tp, $) = (cp', $'), defined 
as follows. If the particle hits Q at a point 7(y?) under an angle of incidence $ 6 (0, n) with 
the tangent vector at 7(y?), then, as the motion is free inside Q, the next impact point is 7(<^')' 
the intersection point with the boundary and the next angle of incidence is G (0, n), as in 
FigureQ] A straightforward computation shows that f(ip, $) = (ip', if and only if 

|7'(y?)| cos$ = -dih((p,Lp'), W(<p')\ cosi?' = d 2 h(ip, ip'), (6) 

where h : T 2 \ {ip' 7^ (/?}—)• R is given by h(cp,cp') = l'j(p) — 7(v 3 ')l- Besides, the twist 
condition holds: dtp'/dfi = h((p, <p')/\ 7' (<p>') \ sin •&' > 0. Finally, it is geometrically clear that 
/ verifies the rigid boundary conditions with w_ = and co + = 2ix. 

A remark is in order. Equations © differ slightly from equations ©, but identity © still 
holds and so the theory developed in the previous section still applies. 
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Figure 2. Left: A (1, 4) -resonant convex smooth caustic C. Right: Its two smooth RICs 

T~ = graph i? - and Y~ = graph i)~ in the phase space Z = T x (0, 7r). 

Obviously, one could write the map in the canonical coordinates — arclength parameter 
for the boundary and cos $ as its conjugate — in order to have h as a generating function, but 
this is not a wise choice when dealing with ellipses. 

Let us assume that there exists a closed convex smooth caustic C contained in the region 
enclosed by Q. Then the billiard map / : Z — > Z has two smooth RICs T 1 * 1 = graph ^ C Z. 
The functions 'i9 ± : T — ?• (0, 7r) are easy to understand: $ + (y?) and #"(<£>) are the angles 
determined by the two tangent lines to the caustic C from the point 7(y?) G Q, see Figure[2] In 
particular, $~ ((p)+ r d + (Lp) = it. To fix ideas, we will assume that Y~ and T + correspond to the 
billiard motion around C in the couterclockwise and clockwise senses, respectively. Hence, 
< $~{(p) < vr/2 < $ + (y?) < n. There is an explicit formula relating the parametrization 
of the billiard curve Q, the parametrization of the caustic C, and the functions ■& ± , See, for 
instance, IfTlfTOl. 

Let Q be a closed strictly convex smooth billiard boundary with a (m, n) -resonant convex 
caustic C, so that its RIC Y~ is (m, n)-resonant and its RIC T + is [n — m, n)-resonant. We 
say that C persists under a perturbation Q e = Q + 0(e) whenever the perturbed billiard curve 
has a (m, n)-resonant caustic C e for any small enough e such that C £ = C + 0(e). 

Let f e be the billiard map inside Q e and Li(ip) and L^(lp) be the subharmonic Melnikov 
potentials of the resonant RICs Y~ and T + under the area-preserving twist perturbation f e . 
Both potentials coincide, due to the time reversibility of the billiard dynamics. Therefore, we 
can skip the ± signs. In this context, we will say that Li(ip) is the subharmonic Melnikov 
potential of the resonant caustic C for the perturbation Q t . 

Corollary 9. IfLi(<p) is not constant, then the resonant caustic C does not persist under the 
perturbation Q e . 
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From now on, we will assume that the unperturbed billiard boundary is the ellipse 



,2.^ 2 , y 2 



Q= -{q=(x,y) eR 2 : — + »— = l\, a>b>0 



a 



It is known that the convex caustics of the billiard inside Q are the confocal ellipses 

Cx = {g = (x, y) e R 2 : + = l} , < A < b. 

Let p(\) be the rotation number of the elliptical caustic C\. Then p : (0,6) — » E is an 
analytic increasing function such that p(0) = and p(b) = 1/2. See, for instance, Hi. 
Thus, there is a unique (m, n)-resonant elliptical caustic for any relatively prime integers m 
and n such that 1 < m < n/2. We shall see that the caustic parameter A G (0,6) of the 
(m, n)-resonant caustic is implicitly determined by means of an equation containing a couple 
of elliptic integrals, see equation (flOl) . 

The following lemma on elliptic billiards is useful to simplify the expression of the 
subharmonic Melnikov potential later on. 

Lemma 10. Let (qj)j^z be any billiard trajectory inside the ellipse Q with caustic C\. Let 
Pj = (qj+i — qj)/\qj + i — qj \ be the unit inward velocities of the trajectory. Then 

ablgj-x - pj, D~ 2 qf) = 2A, Vj G Z, 

where D = diag(a, b) is the diagonal matrix such that Q = {q G M 2 : (q, D~ 2 q) = 1}. 

Proof. We shall prove that given any point q = (x,y) G Q and any unit inward vector 
p = (u, v) G S 1 , the line i = {q + rp : r G M} is tangent to the conic C\ if and only if 

A = —{bxu/a + ayv/b) = —ab(p, D~ 2 q). 

To begin with, we note that the line £ is tangent to the conic C\ if and only if the equation 
of second order in the variable r given by 

(x + ru) 2 /(a 2 - A 2 ) + (y + rv) 2 /(b 2 - A 2 ) - 1 = 

has zero discriminant, which is equivalent to the equation 

xu yv A 2 / u 2 v 2 \ / x 2 y 2 

1 O 1 o I O \ O TO AO J I O AO 1 O AO 1 



a 2 - A 2 b 2 - X 2 J \a 2 - A 2 b 2 - A 2 J \a 2 - A 2 b 2 - A 2 
After some simplifications, we can rewrite this equation as 

(xv - yu) 2 = {b 2 - A V + (a 2 - A> 2 = aV + b 2 u 2 - A 2 , 

since u 2 + v 2 = 1. Next, using that x 2 /a 2 + y 2 /b 2 = 1, we obtain that 

A 2 = (aV + b 2 u 2 )(x 2 /a 2 + y 2 /b 2 ) - (xv - yu) 2 = (bxu/a + ayv/b) 2 . 

Thus, we have two possibilities: A = ab(p, D~ 2 q) or A = —ab(p, D~ 2 q) . The first one is 
discarded, because A > and (p, D~ 2 q) < 0. The second inequality follows from the fact 
that the vector p points inward Q at q, whereas D~ 2 q is an outward normal vector to Q at q. 
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Finally, we note that ~Vj-i = (<?j-i - <?j)/l<?j-i ~ Qj\ andp 3 - = (q j+ i - qj)/\q j+ i - qj\ 
are the two unit vectors that point inward Q at the impact point g, and give the two tangent 
directions to the caustic C\. Therefore, A = o6(pj_i, D~ 2 qf) = —ab{pj, D~ 2 qf). □ 

Proposition 11. Let C\ be the (m,n)- re sonant elliptical caustic confocal to the ellipse Q. 
Given any angle ip G T, let qj = (a cos ipj, b sin <pj) be the vertexes of the (m, n)-gon inscribed 
in Q and circumscribed around C\ such that q = (a cos (p, b sin ip). Then the subharmonic 
Melnikov potential of the caustic C\for the perturbed ellipse ([/]) is 

n-l 

L 1 (<p) = 2\J2^M- (V) 

j=0 

Proof. The parametrization of the perturbed ellipse CQ) is given by 

%(<p) = (ccosh/i e (v?) cos<y9, csinh/i e (v?) simp) = 7o(v?) + Hiif) + 0( e2 )> 

where 7o(y?) = (a cos <p, 6 sin ip), ji(ip) = ab[ii(<p)D~ 2 j ((p), and D = diag(a, b) as above. 
The generating function of the billiard map inside the perturbed ellipse is 

h £ ((p,(p') = \j e ((p') - 7e(^)l = h Q ((p,(p') + ehi(<p,<{/) + 0(e 2 ). 

The first terms of this expansion verify the identities ho(ip,ip') = |7o(v 9 ') — 7o(v)| an d 
h ((p,(p')hi((p,(p') = (7o(^) -7o(^),7i(<^') ~7i(y?))- 

Let (qj)jez De me billiard trajectory inside the ellipse Q with caustic C\ such that 
qj = 7o(yj) and ip = (p. The unit inward velocities of this trajectory are 

Qj+i ~ Qj = ToQgj+i) - 7o (<Pj) 
3 h {ipj,(p j+ i) 
It follows from Proposition [7] that the subharmonic Melnikov potential is 

n-l 

3=0 
n-l 

= 5^(Pj,7i(^j+i) -7i(Vj)> 

j=0 
n-l 

= ^^{p^/ii^+O-D^gi+i - m{ipj)D- 2 qj) 

3=0 

n-l 

3=0 

n-l 

= 2A^/ii(^). 

i=o 

We have used the periodicity in the fourth equality and Lemma [101 in the last one. □ 
Next, we give a couple of sufficient conditions for the subharmonic Melnikov potential 
to be constant. These conditions are trivial. Nevertheless, they play a key role in our problem. 
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Concretely, we shall check later on that they are also necessary conditions in the class of 
2/T -periodic entire functions //i (<£>). 

Corollary 12. Let fii((p) be any 2n-periodic smooth function. 

(i) If the period n is odd, then fii(<p) constant =>■ L\(ip) constant. 

(ii) If the period n is even, then ^[(f) n-antiperiodic =>• L\{ip) constant. 

Proof. The case n odd is obvious. If n is even, the (m, n)-gons inscribed in Q and 
circumscribed around C\ are symmetric with respect to the origin, so <Pj+ n /2 = fj + tt and 

n-l re/2-1 

L' 1 (^) = 2A^ / i' 1 (^) = 2A ]T (^O+^+t))- 
j'=o j=o 

In particular, n even and ^[((p) 7r-antiperiodic =>- -^(v 9 ) = -^(v 9 ) constant. □ 
The subharmonic Melnikov potential of the (m, n) -resonant caustic for the perturbed 
circle © is 

n-l 

LxiB) = 2r sm(m7c/n) ^ n(^), % = + 27rmj/n, (8) 

i=o 

see [16, Proposition 10]. We recall that A = r sin(m7r/n) is the (m, n) -resonant caustic 
parameter of the circle of radius r . Besides, all the (m,n)-gons inscribed in the circle of 
radius r and circumscribed around the circle of radius A = r sin(m7r/n) are regular, so 
their vertexes are of the form qj = (r cos 9 j,r sin 9 j) with 9j = 9 + 2ixmj/n. Hence, the 
function © is the limit of function CO) when both a and b tend to r . 

Although functions © and ([8]) look quite similar, they hide a crucial difference. There 
is a simple formula for the 9j angles, but not for the cpj ones. This has to do with the 
fact that the billiard trajectories inside a circle of radius r sharing a circular caustic with 
radius A = r sm(<5/2) have a rigid angular dynamics of the form 9 i-> 9 + 5. On the 
contrary, such a rigid angular dynamics does not take place for elliptic tables when the angle 
cp is considered, which is a source of technical difficulties in the study of the subharmonic 
Melnikov potential ©. Nevertheless, it is possible to define a new angular parameter t over 
the ellipse Q in such a way that all billiard trajectories inside Q sharing the elliptical caustic 
C\ have a rigid angular dynamics of the form t h-> t + 5, for some constant shift 5 = 5(A). 

We need some notations on elliptic functions in order to define this angular parameter t. 
We refer to [ffllHl for a general background on elliptic functions. Given a quantity k 6 (0, 1), 
called the modulus, then K = K{k) = J^ 2 (l — k 2 sin 2 0)~ 1/,2 d0 is the complete elliptic 
integral of the first kind. We also write K' = K'{k) = K(yl — k 2 ). The amplitude function 
Lp = am t is defined through the inversion of the integral 

t= f V (1 - k 2 sin 2 (f))- 1/2 d(j). 
Jo 

Then the elliptic sinus and the elliptic cosinus are defined by the trigonometric relations 

snt = siny), cnt = cos ip, 
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respectively. Dependence on the modulus is denoted by a comma preceding it, so we can write 
am(t, k), sn(i, k), and cn(t, k) to avoid any confusion. In the following lemma it is stated that 
the angular dynamics becomes rigid in the angular parameter t given by tp = am(t, k). It 
suffices to find the suitable modulus k for each elliptical caustic C\. 

Lemma 13. Once fixed any caustic parameter A G (0, b), we set the modulus k G (0, 1) and 
the constant shift 5 G (0, 2K) by the formulae 



where d G (0, n) is the angle such that sin($/2) = A/6. Let 

qj = (a cos ipj, frsin <pj) = (acn(tj,k),bsn(tj,k)) 
be any billiard trajectory inside the ellipse Q with caustic C\. Then tj + i = tj + 5. 

Proof. By definition, cpj = &m(tj,k), so tj + \ — tj = J^ J+1 (1 — k 2 sin 2 0) _1 / 2 d0. These 
integrals are equal to a constant 5 that depends only on C\, see B51 page 1543]). The formula 
for the constant shift is given in [5, page 1540]. □ 

Remark that if a = b = r then the modulus k is equal to zero, the complete elliptic 
integral K is equal ton/2, the amplitude function is the identity, the elliptic sinus/cosinus are 
the usual sinus/cosinus, the shift 5 G (0, n) is given by A = r sin(5/2), and the dynamical 
relation — tj + 5 becomes (pj+i = <fj + Thus, we recover the known rigid angular 
dynamics for circular tables as a limit of the formulae for elliptic tables. 

From now on, k and S will denote the modulus and the constant shift defined in ©. Thus, 
we shall skip the dependence of the elliptic functions on the modulus. We note that C\ has 
eccentricity k. Besides, C\ is the (m, n)-resonant elliptical caustic if and only if 



This identity has the following geometric interpretation. When a billiard trajectory makes 
one turn around C\, the old angular variable ip changes by 2n, so the new angular variable 
t changes by AK. On the other hand, we have seen that the variable t changes by 5 when a 
billiard trajectory bounces once. Hence, a billiard trajectory inscribed in Q and circumscribed 
around C\ makes exactly m turns around C\ after n bounces if and only if (flOl) holds. 

Proposition 14. Let fJ>i(<f) be any 2n -periodic entire function. 

(i) If the period n is odd, then L\{ip) constant ^ P-i(^>) constant. 

(ii) If the period n is even, then Li(tp) constant <=>■ fi[((p) n-antiperiodic. 

Proof. Let A = 2K + 2K'i and z(t) = cn t + i sn t. If ip = am t, then 





(9) 



nS = AKm. 



(10) 



c 



cos (p + isiny? = cnt + isnt = z(t), 
cos ip — i sin (p = cn t — i sn t = z(t + A) . 
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We have used that the elliptic cosinus is A-periodic, but the elliptic sinus is A-antiperiodic. 
We also recall that the elliptic cosinus/sinus are 2i^-antiperiodic meromorphic functions on 
the whole complex plane whose unique singularities are the points of the form 

r rjS = 2Kr + (1 + 2s) K'i, r, s G Z. 

Besides, these singularities are just simple poles whose residues are 

res(cn;r riS ) = (-l) r+s+1 i//c, res(sn; r r , s ) = (-l) r /k. 

Thus, z(t) is a 2i^-antiperiodic meromorphic function whose unique singularities are the 
points of the set 

P = {r r ,2s+i :r,seZ} = n + 2KZ + AK'iZ, t, = r 0j _i = -K'i. 

As before, these singularities are just simple poles. 

Let J2iez P>i elllf be the Fourier expansion of //i (<£>). Then 

yUi(amt) = nx{ip) = 22flie ll(p = p,-{z(t + A)) + fi + fi + (z(t)), 

where p,-(z) = Xlz^i A-z- 2 ' an d A+( z ) = X^i/W- We note that the functions £l±(z) are 
entire, because fj,i((p) is entire. Besides, 

n-1 

Li(amt) =X x (yj) = 2A^/i 1 (^ i ) = 2A (L_(t) + n/x + £+(*)) > ( n ) 

i=o 

where = X}j=o ft-( z (t + A + and = YllZo f^+( z (t + jS)). Let us study the 

behaviour of these two functions around the point t* = —K'i. Concretely, we shall prove that 
L-(t) is analytic at t = r*, whereas L+(t) has a nonremovable singularity at t = r* provided 
is nonconstant and n is odd, or provided /^(v?) is not 7r-antiperiodic and n is even. 
We begin with a couple of simple observations. If j G {0, . . . , n — 1}, then: 

a) 3(r* + A + jS) = K', so r* + A + jS & P; and 

b) n +j5ePo AKmj/n = j5 G 2KZ <=> 2jm G nZ <=> 2j GnZ^jG {0, n/2}. 
Here, we have used that 5 G R, equation (flOl) . and gcd(m,n) = 1. Besides, we stress 
that the equality j = n/2 only can take place when n is even. 

We deduce the following results from the above observations. 

1) L-(t) is analytic at t = r*, because so are z(t + A + jS) for j = 0, . . . , n — 1. 

2) If n is odd and (</?) is nonconstant, then: 

- The function (l+(z) is nonconstant and entire; 

- The function L + (t) — fi + (z{t)) = Y^jZx fa+( z (t + jS)) is analytic at t = r*; 

- The composition fi + (z(t)) has a nonremovable singularity at t = r*; and 

- The function (fTTT) is nonconstant, since it has a nonremovable singularity at t — r*. 

3) If n is even and ^[(ip) is not 7r-antiperiodic, then: 

- The sum a(z) = fi+(z) + jl + (—z) = 2 Y^tli faiz 21 is a nonconstant entire function; 

- z(t + n5/2) = z(t + 2Km) = (-l) m z(t) = -z{t), since m is odd; 
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- £+(*(*)) + fr+(z(t + nS/2)) = a(z(t)); 

- The function L + (t) — a(z(t)) is analytic at t = t*; 

- The composition a(z(t)) has a nonremovable singularity at t = t*; and 

- The function (flTT) is nonconstant, since it has a nonremovable singularity at t = r*. 

Therefore, the proof follows by combining the above results with Corollary [12] □ 
Finally, we note that our main result (namely, Theorem Q] stated in the introduction) 
follows directly from Corollary [9] and Proposition [141 
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